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^ , Abstract 

The Minimal Supersymmetric extension of the Standard Model (MSSM) can pro- 
vide the correct neutralino relic abundance and baryon number asymmetry of the 
universe. Both may be efficiently generated in the presence of CP violating phases, 
light charginos and neutralinos, and a light top squark. Due to the coannihilation of 
the neutralino with the light stop, we find a large region of parameter space in which 
the neutralino relic density is consistent with WMAP and SDSS data. We perform a 
detailed study of the additional constraints induced when CP violating phases, consis- 
tent with the ones required for baryogenesis, are included. We explore the possible tests 
of this scenario from present and future electron Electric Dipole Moment (EDM) mea- 
surements, direct neutralino detection experiments, collider searches and the 6 — > 57 
decay rate. We find that the EDM constraints are quite severe and that electron EDM 
experiments, together with stop searches at the Tevatron and Higgs searches at the 
LHC, will provide a definite test of our scenario of electroweak baryogenesis in the 
next few years. 



1 Introduction 



The nature of the dark matter and the source of the baryon-anti-baryon asymmetry are two 
of the most important questions at the interface of particle physics and cosmology. Recent 
improvements in the astrophysical and cosmological data, most notably due to the Wilkinson 
Microwave Anisotropy Probe (WMAP) [1] and the Sloan Digital Sky Survey (SDSS) [2], have 
determined the matter and baryon densities of the Universe to be VLm^? = 0.135^q^qI and 

^]^/i2 = 0.0224 ± 0.0009, (1) 

respectively, with h = 0.71^q q3. Together these imply a (dominantly) cold dark matter 
density of 

ncDAih' = 0.11261°:°}^}, (2) 

at 95% CL. Such precise determinations of Qsh"^ and QcDAih"^ impose severe constraints on 
any particle physics model that tries to explain one or both of these values. 

The Standard Model of particle physics (SM) has been tested extensively by collider 
experiments, and has so far withstood all of them. However, the SM performs considerably 
worse when it comes to cosmology, and can account for neither the baryon asymmetry, nor 
the dark matter. Furthermore, in the SM, the electroweak scale is unstable under quantum 
corrections suggesting that an extension of the SM description is required at energies near 
the TeV scale. A particularly attractive way to stabilize the weak scale is to introduce 
supersymmetry [3]. Remarkably, the minimal supersymmetric extension of the standard 
model, the MSSM, can also explain the baryon asymmetry, and contains an excellent dark 
matter candidate in the lightest supersymmetric particle (LSP). 

The LSP of the MSSM is stable if R-parity is imposed. If, in addition, the LSP is 
neutral under SU{3)c x U{1)em, it is a candidate for cold dark matter. One such particle 
is the lightest neutralino. This particle tends to have a mass of order of the weak scale 
and electroweak strength couplings, and therefore naturally gives rise to a dark matter relic 
density close to the measured value. The fact that a stable particle with electroweak strength 
couplings and mass of order 1 TeV naturally generates a relic density near the required value 
can be taken as further motivation for new physics at the TeV scale. 

In general, any mechanism for baryogenesis must fulfill the three Sakharov require- 
ments [1]; namely baryon number (B) violation, CP violation, and a departure from equi- 
librium (unless CPT is violated, see for instance ^). All three requirements are satisfied in 
both the SM and the MSSM during the electroweak phase transition, and this is the basis for 
electroweak baryogenesis (EWBG) [6]. However, as we will discuss below, while electroweak 
baryogenesis may be realized in the MSSM, SM processes cannot generate a large enough 
baryon asymmetry during the electroweak phase transition. 

Baryon number violation occurs in the SM and the MSSM due to anomalous sphaleron 
transitions that violate (B + L) 0. These transitions are exponentially suppressed at low 
temperatures in the electroweak broken phase [8], but become active at high temperatures 
when the electroweak symmetry is restored [9]- In the absence of other charge asymmetries. 
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like (B—L), they produce baryons and anti-baryons such that the net baryon number relaxes 
to zero, and so do not by themselves generate a baryon asymmetry [6l [10] 

If the electroweak phase transition is first order, bubbles of broken phase nucleate within 
the symmetric phase as the Universe cools below the critical temperature. These provide 
the necessary departure from equilibrium. EWBG then proceeds as follows [H]. CP vi- 
olating interactions in the bubble walls generate chiral charge asymmetries which diffuse 
into the symmetric phase in front of the walls. There, sphaleron transitions, which are ac- 
tive in the symmetric phase, convert these asymmetries into a net baryon number. This 
baryon number then diffuses into the bubbles where the electroweak symmetry is broken. 
Sphaleron transitions within the broken phase tend to destroy the baryon number generated 
outside the bubble. To avoid this, the sphaleron transitions within the broken phase must 
be strongly suppressed. This is the case provided the electroweak phase transition is strongly 
first order [12], 

v{T,)/T, > 1 , (3) 

where v{Tc) denotes the Higgs vacuum expectation value at the critical temperature Tc. 

The strength of the electroweak phase transition may be determined by studying the 
finite temperature effective Higgs boson potential. The Higgs vacuum expectation value at 
the critical temperature is inversely proportional to the Higgs quartic coupling, related to 
the Higgs mass. For sufficiently light Higgs bosons, a first-order phase transition can be 
induced by the loop effects of light bosonic particles, with masses of order the weak scale 
and large couplings to the Higgs fields. The only such particles in the SM are the gauge 
bosons, and their couplings are not strong enough to induce a first-order phase transition 
for a Higgs mass above the LEP II bound [IB]. 

Within the MSSM, there are additional bosonic degrees of freedom which can make the 
phase transition more strongly first-order. The most important contribution comes from 
a light stop, which interacts with the Higgs field with a coupling equal to the top-quark 
Yukawa. In addition, a light stop has six degrees of freedom, three of colour and two of 
charge, which further enhances the effect on the Higgs potential. Detailed calculations show 
that for the mechanism of electroweak baryogenesis to work, the lightest stop mass must be 
less than the top mass but greater than about 120 GeV to avoid colour-breaking minima. 
Simultaneously, the Higgs boson involved in breaking the electroweak symmetry must be 
lighter than 120 GeV [ll]-[2T], only slightly above the present experimental bound [22], 

rriH > 114 GeV, (4) 

which is valid for a Higgs boson with SM-like couplings to the gauge bosonsQ 

The combined requirements of a first-order electroweak phase transition, strong enough 
for EWBG, and a Higgs boson mass above the experimental limit severely restrict the allowed 
values of the stop parameters. To avoid generating too large a contribution to Ap, the light 
stop must be mostly right-handed. Since the stops generate the most important radiative 

^The requirements of a light stop and a light Higgs boson may be relaxed in non- minimal supersymmetric 
extensions. See, for instance, Refs. [23]-[30]. 
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contribution to the Higgs boson mass in the MSSM [31], the other stop must be considerably 
heavier in order to raise the Higgs mass above the experimental bound, Eq. (jl]). For the stop 
soft supersymmetry-breaking masses, this implies [18] 

ml;^ < 0, (5) 
^Q3 ^ (ITeVf. 

A similar tension exists for the combination of soft SUSY breaking parameters defining the 
stop mixing, \At—^*/ tan P\/mQ^, and tan p. Large values of these quantities tend to increase 
the Higgs mass at the expense of weakening the phase transition or the amount of baryon 
number produced. The allowed ranges have been found to be [IB] 

5 < tan/? < 10, (6) 
0.3 < \At-fi*/tan(3\/mQ, < 0.5. 

A strong electroweak phase transition is only a necessary condition for successful EWBG. 
In addition, a CP violating source is needed to generate a chiral charge asymmetry in the 
bubble walls. Within the MSSM, the dominant source is produced by the charginos, and 
is proportional to /m(/iM2) [321 133]. For this source to be significant, the charginos must 
be abundant in the plasma, which requires that they not be too much heavier than the 
temperature of the plasma, T ~ T^. In the recent analysis of Ref. [33], the authors found 
the bounds 

\ArgiixM2)\ > 0.1, (7) 
/i. Ma < 500 GeV . 

These conditions are very relevant to the issue of neutralino dark matter. 

The need for a large CP violating phase, Eq. ([7]), implies that there is a danger of 
violating the experimental bounds on the electric dipole moments (EDM) of the electron, 
neutron, and ^^^Hg atom since phases generate new contributions to these EDM's. The 
leading contributions arise at one loop order, and they all contain an intermediate first 
or second generation sfermion. They become negligible if these sfermions are very heavy, 

> 10 TeV. Such large masses have only a very small effect on EWBG. At two-loop order, 
if Arg{n M2) 7^ 0, there is a contribution involving an intermediate chargino and Higgs 
boson [311 [35]. Since EWBG requires that this phase be non-zero and that the charginos be 
fairly light, the two-loop contribution is unavoidable if EWBG is to be successful. Thus, EDM 
limits strongly constrain the EWBG mechanism in the MSSM. Similarly, the branching ratio 
for b —>■ s'-f decays is also sensitive to this phase, and therefore imposes a further constraint 
on the EWBG mechanism. 

In a previous work [36j, some of the present authors investigated the neutralino relic 
density in the presence of a light squark, as required for EWBG, but without including the 
effects of CP violating phases in the calculations. Here, we extend the analysis to study in 
detail the effect of phases in order to better understand the relationship between EWBG 
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and dark matter within the MSSM. The outhne of the paper is as follows. In Section [2] we 
investigate the relic density of a neutralino LSP in the presence of both a light stop and 
CP violating phases. Section [3] examines the prospects for direct detection of the neutralino 
dark matter in laboratory experiments, again including CP violating phases. In Section HI 
we will look at the constraints on the phases needed for EWBG due to the electron EDM 
and flavour- violating b s'-f transitions. Finally, Section [5] is reserved for our conclusions. 

2 Neutralino Dark Matter 

As discussed in the introduction, the dual requirements of successful EWBG and a lightest 
Higgs boson with mass greater than the LEP II bounds strongly constrain the parameter 
space of the MSSM [37]. One of the stops must be light, with mass less than that of the 
top, and mostly right-handed. Furthermore, the charginos must not be too heavy, and 
the combination yuM2 must have a non-negligible phase. These conditions have important 
implications for neutralino dark matter. 

First of all, if the lightest neutralino is to be the source of the observed dark matter, it 
must be lighter than the light stop so that it be stable. Secondly, in much of the parameter 
space of interest the light stop is only slightly heavier than the neutralino LSP implying that 
stop-neutralino coannihilation is significant. Finally, a phase for fiM2 modifies the masses of 
the neutralinos and their couplings to other particles, and can also affect the relative phase 
between the various contributions to the annihilation cross-section. The effect of CP violating 
phases on neutralino dark matter has been considered previously by several groups [40j- 
|43] . However, in all of these analyses the regions of MSSM parameter space considered were 
much different from the restricted subset required for EWBG, and in particular, none of 
them included a light stop. 

To simplify the analysis, we shall assume throughout this work that the gaugino mass 
parameters Mi and M2 are related by the standard unification relation, M2 = {sll 9i) — 
2 Ml. The stop soft parameters are largely fixed by the EWBG and Higgs mass conditions. 
We take them to be 

ml^ ^ 

niQ, = 1.5 TeV (8) 
\Xt\ = \At- jj*/ tan (3\ = 0.7 TeV. 

We also set m^g = tjil^ = rriE^ = 1 TeV. EWBG and the Higgs mass constraint also require 
5 < tan /3 < 10 and Ma > 200 GeV. For concreteness, we shall consider the values 

tan/3 = 7 (9) 
Ma = 200, 1000 GeV. 

The first and second generation sfermion soft masses are taken to be very large, mj >10 TeV. 
As we will discuss in section HJ this is necessary to avoid the electron, neutron, and ^^^Hg 
EDM constraints in the presence of large phases. 
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The only phase that we consider in this work is the one directly related to EWBG, namely 
Arg{n M2). We will assume further that this phase is the result of a common phase for the 
gaugino mass parameters. With this assumption, all CP violating effects are confined to 
the chargino and neutralino sectors, or the loop corrections induced by them|l By means of 
a U(l)ii; transformation, we may transfer the gaugino phase into the fi parameter and the 
trilinear Af terms. Under this transformation, the effective values of these parameters are 
shifted according to 

M, ^ M^e-^^ (10) 
Af Afe-^\ 

with the remaining MSSM parameters left unchanged. For consistency of notation with [33], 
we will implicitly make a U(1)r rotation such that the gaugino masses are all real and 
positive, and the [i parameter and the A^ terms have equal and opposite phases (up to a 
possible relative sign). 

As a further simplification, we will neglect the mixing between CP-even and CP-odd 
Higgs bosons due to these phases. While this mixing can be significant in some regions of the 
MSSM parameter space, especially for large values of tan /3, |yu| and \A^, and small M^^ [55] , 
we have checked that the mixing (induced by chargino and neutralino loop corrections) is 
small (< 3%) for the parameters considered here, where tan/5 takes only moderate values 
and the only relevant phase is the one associated with the gaugino sector. We also note that 
in [l2] the effect of Higgs mixing on the neutralino relic density was found to be small, even in 
the large tan/5 regime, where the Higgs boson mixing is much larger. The supersymmetric 
corrections to the bottom mass [39] are also suppressed in the region of parameter space 
considered here, and hence all relevant CP violating effects are associated with the tree-level 
effect on the neutralino masses and couplings. 



2.1 Relic Density 

We compute the relic abundance of neutralinos by numerically solving the Boltzmann equa- 
tion, 

for the number density of the supersymmetric particles n = YliLi ^i- -Due to conservation of 
R parity, the present value of n is equal to the number density of the lightest neutralino ni. 
In Eq. (fTTj) H = lOO/i km/sec/Mpc, neq is the value of n at thermal equilibrium, and 

WeffV){x) = '——^ (12) 



4 / Y^A^ 

X 



^We do not consider the effects of a gluino phase. For the parameters considered in the present work, we 
expect that such a phase would only have a very small effect. 
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is the thermally averaged annihilation cross section. This quantity is a function of x = rrii/T, 
and is given in terms of the individual annihilation cross sections o"jj(a) of the processes ij 
SM and/or Higgs particles. The energy and mass fractions a = ^/s/ml and bi = mj/mi also 
enter via X{a'^, bf, b]) = a'^ + bj + b'^ - 2{a%j + a^b] + bjb]). In Eq.([l2l) gi is the number of 
degrees of freedom of the iih supersymmetric partner, and Ki is the modified Bessel function 
of the second kind of order I. The mass of the lightest neutralino is denoted by mi. 

In our calculation all relevant annihilation and coannihilation processes are included as 
described in Ref . |44j . Besides neutralino self-annihilations, coannihilations of the lightest 
neutralino with the lightest stop and the lighter chargino, and annihilations of the lightest 
stop and chargino effect significantly our numerical results. The complex phases enter our 
relic density calculation directly through the couplings and indirectly through the masses of 
the neutralinos and charginos. After diagonalization of the gaugino and sfermion complex 
mass matrices, we calculate the annihilation cross sections with complex couplings. In doing 
this, we follow techniques used in Refs. [HI HB] . 




a^, = 300 30 3xl0''°pb □ Q.I1 < 0.095 ^ m^,< 103.5 GeV a„. = 300 30 3 xlO''°pb 



r I 0.095 <Q.h^ < 0.129 

Figure 1: Neutralino rehc density for Ma = 200 GeV (left) and Ma = 1000 GeV (right), 
and Arg{fi) = 0. 

Figures [IH show the dependence of the neutralino relic density on |yu| and Mi for tan (5 = 
7, Ma = 200 GeV (left) and Ma = 1000 GeV (right), and three values of the /x phase: 
Argyll) = 0, 7r/2, tt. Values of the phase equal to or vr are representative of what happens 
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for small phases, like the ones consistent with the generation of the baryon asymmetry when 
~ M2 and M4 < 300 GeV where there is a resonance in the amount of baryon number 
produced [33]. On the other hand, large values of the phase, close to 7r/2, tend to be 
necessary to generate the baryon asymmetry outside of the resonant region, particularly for 
large values of Ma, for which the EDM constraints become less severe. 

The green (medium gray) bands in Figures [TH show the region of parameter space where 
the neutralino relic density is consistent with the 95% CL limits set by WMAP data. The 
regions in which the relic density is above the experimental bound and excluded by more 
than two standard deviations are indicated by the red (dark gray) areas. The yellow (light 
gray) areas show the regions of parameter space in which the neutralino relic density is less 
than the WMAP value. An additional source of dark matter, unrelated to the neutralino 
relic density, would be needed in these regions. Finally, in the (medium-light) gray region 
at the upper right the lightest stop becomes the LSP, while in the hatched area at the lower 
left corner the mass of the lightest chargino is lower than is allowed by LEP data [17]. 




a^, = 300 30 3xl0''°pb □ Q.I1 < 0.095 ^ m^,< 103.5 GeV a„ = 300 30 3 xlO''°pb 



= 7_ 10 13xl0'^^ecm □ 0.095 <Q.h^ < 0.129 = 7 1.2 1.4xl0'"ecm 

Figure 2: Neutralino rehc density for Ma = 200 GeV (left) and Ma = 1000 GeV (right), 
and Arg{fi) = 7r/2. 

These figures are qualitatively similar, but do show some differences due to the change in 
the phase of fi. Before discussing the effect of the phase, we will examine the general features 
of Figures [T]l3l For Ma = 1000 GeV and for all three phase values, the region where the relic 
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density is too high consists of a wide band in which the hghtest neutrahno has mass between 
about 60 and 105 GeV and is predominantly Bino. Above this band, the mass difference 
between the neutrahno LSP and the hght stop is less than about 20 GeV, and stop-neutralino 
coannihilation as well as stop-stop annihilation are very efficient at reducing the neutrahno 
abundance. For Ma = 200 GeV, instead, the contribution to neutrahno annihilation from 
s-channel exchange of heavy CP-even and CP-odd Higgs bosons is enhanced by a resonance 
around m^^ ^ 100 GeV. This restricts the band in which the relic density is too high to 
the region where the lightest neutrahno has mass between about 60 and 85 GeV, and is also 
mostly Bino. For both values of Ma-, there is an area below the disallowed band in which the 
neutrahno mass lies in the range 40-60 GeV, and the neutrahno annihilation cross-section is 
enhanced by resonances from s-channel and exchanges. 




CT^,, = 300 30 3xl0''°pb □ Q.h^< 0.095 ^ m^,< 103.5 GeV a„ = 300 30 3 xlO''°pb 



H 0.095 <£lh^ < 0.1 29 

Figure 3: Neutrahno rehc density for Ma = 200 GeV (left) and Ma = 1000 GeV (right), 
and Arg{^) = ir. 

The relic density is also quite low for smaller values of In these regions, the neutrahno 
LSP acquires a significant Higgsino component allowing it to couple more strongly to the 
Higgs bosons and the Z^. For Ma = 1000 GeV, this is particularly important in the region 
near (|/i|,Mi) = (175,110) GeV where the neutrahno mass becomes large enough that 
annihilation into pairs of gauge bosons through s-channel Higgs and Z° exchange and t- 
channel neutrahno and chargino exchange is allowed, and is the reason for the dip in the 
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relic density near this point. Since the corresponding couphngs to the gauge bosons depend 
on the Higgsino content of the neutrahno, these decay channels turn off as |/i| increases. For 
higher Mi values, the lightest neutralino and chargino masses are also close enough that 
chargino-neutralino coannihilation and chargino- chargino annihilation substantially increase 
the effective cross section. 

In Figures ([UlS]), we have taken M2 = {g2/gf)Mi, as suggested by universality. Because of 
this, smaller values of Mi and n are excluded by the lower bound on the chargino mass from 
LEP data [17], as indicated by the hatched regions in the figures. This constraint becomes 
much less severe for larger values of the ratio M2 / Mi . We also find that increasing this ratio 
of gaugino masses (with Mi held fixed) has only a very small effect on the neutralino relic 
density. 



For the parameters considered in the previous section, relevant for EWBG within the MSSM, 
CP violating phases modify the values of the neutralino relic density but have only a mild 
effect on the general qualitative features of the allowed parameter space. This is somewhat 
misleading, however, since the value of the relic density at a given point in the \fi\—Mi plane 
can vary markedly with Arg{n). 

The most important effect of varying Arg{fi) is to shift the mass of the neutralino LSP. 
The dependence of the lightest neutralino mass on this phase is shown in Figure H] for 
tan/3 = 7 and three sample values of (|//|,Mi): (|/i|,Mi) = (350,110) GeV, (300,60) GeV, 
and (175, 110) GeV. For M^ = 1000 GeV, these three points are representative of the 
regions where the annihilation cross section is dominated by stop-neutralino coannihila- 
tion ((350, 110) GeV), Higgs boson s-channel annihilation ((300,60) GeV), and annihilation 
into pairs of gauge bosons ((175, 110) GeV). In all three cases, the neutralino mass increases 
with Arg{fj,), by about 3%, 7%, and 11%, respectively. Such a mass shift can significantly 
modify the relic density at a single point where neutralino annihilation is enhanced by a 
resonance or coannihilation with another species. The effect on the net distribution of relic 
densities, on the other hand, is fairly small; shifting the phase tends to translate this distri- 
bution down and to the left in the |/u|-Mi plane. 

The neutralino-Higgs couplings are also quite sensitive to Arg{fi). The couplings of the 
Higgs bosons to a pair of neutralinos are given in [71] , and have the form 



where Pl,r = (l=F75)/2 are the usual chiral projectors. Using these vertices, the spin-summed 
and squared matrix elements for ZiZi ff annihilation via s-channel Higgs exchange are 
proportional to 



2.2 Effects of CP Violating Phases 



ZiZi h^H' 
ZiZi A 



i{FPl + F*Pr) 
i{GPl-G*Pr) 



(13) 




AO 



(14) 
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Figure 4: Mass of the lightest neutrahno as a function Arg{fi) for tan (3 = 1 and three sample 
values of Mi). 

In calculating the thermal average, one integrates these matrix elements over s through the 
range [4m^,cxD) with a Boltzmann factor, Eq. (|T2l) . The Boltzmann suppression is strong 
for a cold relic, so the integral is dominated by the region s ~ 4m^. In particular, this 
means that the terms in Eq. ([T3j) proportional to s have the potential to give a much larger 
contribution to the thermal average than those proportional to (s — 4m^). 
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Figure 5: Variation of the real and imaginary parts of the ZiZih^ coupling with Arg{fi) for 
(|/i|,Mi) = (300,60) GeV (solid and dotted), and (|/x|,Mi) = (175,110) GeV (dashed and 
dash-dotted) . 

The dependence of the ZiZi coupling on Arg{fi) for Ma = 1000 GeV, and (|/i|. Mi) = 
(300,60) GeV and (|/i|,Mi) = (175,110) GeV is shown in Figure [51 Both the real and 
imaginary parts of the couplings are larger in the (|/i|,Mi) = (175, 110) GeV case since for 



10 



these values of the parameters, the neutrahno LSP has a much larger Higgsino component 
than for (|/i|,Mi) = (300,60) GeV, when the neutralino is mostly Bino. The couplings for 
(|yLi|,Mi) = (350,110) GeV, where the LSP is also mostly Bino, are very similar to those 
for (|/i|,Mi) = (300,60) GeV. Setting Ma = 200 GeV has only a small effect on these 
couplings. For both points shown in Figure [6], the imaginary part of the coupling vanishes 
when /i is real, and is largest when /i is pure imaginary, Arg{fi) = 7t/2. The real part of the 
coupling also tends to decrease with Arg{fi) due to an accidental cancellation of terms. This 
behavior may be seen by comparing the region Mi < 60 GeV in Figures [T|, [21 and [3l where 
s-channel exchange tends to be dominant. The relic density in this region is lowest when 
Arg{n) = n/2, Figure Ej while in Figure [31 corresponding to Arg^jj) = vr, the contribution 
from h'^ exchange is much smaller than for other values of this phase. 




0.5 I I.S 2 2.5 3 0.5 I I.S 2 2.5 3 



Arg( \l) Arg( \l) 



Figure 6: Variation of the real and imaginary parts of the ZiZiH^ (left) and Z\Z\J^ (right) 
couphngs with Arg{[i) for (|/i|,Mi) = (300,60) GeV (sohd and dotted), and (|/i|,Mi) = 
(175, 110) GeV (dashed and dash-dotted). 

The couplings of the and bosons to neutralinos are shown in Figure [6] for = 
1000 GeV, and (|/i|,Mi) = (300,60) GeV and (|/i|,Mi) = (175,110) GeV. As with the h° 
coupling, these couplings are nearly unchanged when = 200 GeV, and the couplings 
for (|/i|,Mi) = (350,110) GeV are very similar to those for (|/i|,Mi) = (300,60) GeV. The 
imaginary part of the and A^ couplings vanishes for Arg{^) = 0, vr and is largest near 
Arg{fi) = tt/2, while the real parts of these couplings are largest for Arg{fi) = 0, vr and 
nearly zero when Arg{n) = 7r/2. From Eqs. ([T3| UM . this implies that the contribution of 
s-channel exchange to neutralino annihilation is largest when Arg{^) = and smallest 
for Arg{fi) = 0, vr, and that the opposite is true for s-channel A^ exchange. Interestingly, 
the sum of the A^ and contributions is nearly independent of the phase. We expect this 
to be the case whenever M\ ^ M|, and the heavy CP-even and CP-odd Higgs states are 
nearly degenerate. The same effect was found in [12]. 

We have also investigated the phase dependence of the Zitt coupling which generates the 
most important contributions to stop-neutralino coannihilation. While this coupling does 
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vary somewhat with the phase, the effect of the phase on the neutrahno mass is much more 
important. This is because the coannihilation contribution to the rehc density is suppressed 
by a factor e~^^^~'^^^^^'^^ , where Tj ~ /20 is the neutrahno freeze-out temperature, 
making it very sensitive to the neutrahno mass. 

3 Direct Detection of Dark Matter 

If space around us is filled with relic neutralinos, then it is plausible to try to observe them. 
Indeed, the search for weakly interacting massive particles is in progress via detection of their 
scattering off nuclei by measuring the nuclear recoil. Since neutralinos are non-relativistic 
they can be directly detected via the recoiling off a nucleus in elastic scattering. There are 
several existing and future experiments engaged in this search. These include solid state 
germanium, ionization based detectors such as IGEX |18], HDMS gO], CDMS [50], EDEL- 
WEISS pT] and GENIUS [52]. Solid crystal or liquid Nal based scintillator detectors are 
used for example by DAMA [53] and ZEPLIN [SH [SHI (SH |57] . Liquid, gas or hybrid xenon 
based detector is used by experiments as XENON [58] and UKDMC [59]. Gas target projec- 
tion chambers are utilized in DRIFT ^60j, and metastable particle detectors in SIMPLE [61j 
and PICASSO (621. 
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Figure 7: Spin-independent neutralino-proton scattering cross section as the function of 
Arg{fx), for |/i| = 350 GeV and Mi = 110 GeV, and for mA= 200 (1000) GeV for the upper 
(lower) curve. 

The elastic scattering interactions of neutralinos with nuclei can be described by the sum 
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of spin independent {Cgi ) ^^"^ ^P^^ dependent {Cgf)) Lagrangian terms: 

^-elastic ~ *-SI "I" *-SD- \^^) 

For heavy nuclei the spin independent (SI) cross section, being proportional to the squared 
mass of the target nucleus, is highly enhanced compared to the spin dependent one. For 
the case of a target containing the isotope ^^^I, for example, the enhancement factor is more 
than 10^. For this reason the experimental limits on the spin independent neutralino-nucleon 
cross sections are considerably stronger. 
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Figure 8: Same as Figure [H but for |yu| = 300 GeV and Mi = 60 GeV. 



In what follows, we will focus on the spin independent interactions of neutralinos with 
nuclei. At the parton level, these are mediated by t-channel Higgs and s-channel squark 
exchanges. (Here, we only consider the, so called, scalar contribution and neglect the higher 
order tensor contribution originating from loop diagrams.) The differential scattering rate 
of a neutralino off a nucleus with mass mx takes the form ^63] : 



dasi 
d\q\' 



TTV 



:[Zfp+iA-Z)U]^F\Qr), 



where q 



mx+m 



is the three-momentum transfer, Qr 



2m]v ' 



(16) 



and F'^{Qr) is the scalar 



nuclear form factor, v is the velocity of the incident neutralino and fp and are effective 
neutralino couplings to protons and neutrons respectively. The same formalism was used in 
Ref. [63j to calculate neutralino-nucleon cross sections, and the reader is directed there for 
further details. Since modern experiments express their limits in terms of the neutralino- 
proton cross section, we calculate and plot this quantity in this work. 
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To study the dependence of the neutrahno-proton cross section on complex phases of 
various supersymmetric parameters, we select a point in the examined parameter region 
where constraints from EWBG, the electron EDM and WMAP are simultaneously satisfied. 
Specifically, we examine values of Ma = 200, 1000 GeV and the same Higgsino and neutralino 
mass parameters chosen before, namely (|//|,Mi) = (175,110) GeV, (350,110) GeV and 
(300,60) GeV. As emphasized before, for Ma = 1000 GeV these points correspond to regions 
in which the annihilation cross section is dominated by weak processes, coannihilation with 
the light stop, and s-channel Higgs exchange, respectively. 



Legend: 
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« 0.095 <Q.h^ < 0.129 
X m^ = 200GeV: 

= &E-28 lE-27 e cm 
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d^ = 3E-28 9E-28 e cm 



Arg(^) 

Figure 9: Same as Figure [71 but for = 175 GeV and Mi = 110 GeV. 

Figures [3-[9] show the neutralino-proton cross section versus the phase of fi for the selected 
parameter space points. The most striking feature of these plots is that the cross section 
is suppressed for non-vanishing phases and, except for (|/i|,Mi) = (175,110) GeV with 
Ma = 1000 GeV, nearly vanishes for a given value of Arg{fi). This behavior follows from 
the phase dependence of the Higgs-neutralino couplings. In our case, t-channel and 

exchange diagrams generate the most important contributions to the spin-independent 
neutralino-nucleon elastic scattering cross section. (We checked that the only relatively light 
squark, the lightest stop, contributes only at the percent level via its s-channel diagram.) 
Furthermore, these contributions depend only on the real (scalar) part of the Higgs-neutralino 
couplings [Ml |65]; Re{F) in the notation of Eq. f|T3l) . The large suppression of the cross 
section for particular values of Arg{fi) is due to zeroes of Re{F). 

Consider first the Ma = 1000 GeV lines in Figures [3-[9l For these. Ma ^ Mh 3> mh, so 
the contribution of the heavier scalar Higgs is suppressed relative to the lighter state, and 
the neutralino-proton scattering is dominated by t-channel exchange. Comparing the real 
part of the /i°-neutralino coupling for (|/i|,Mi) = (300,60) GeV shown in Figure [5] to the 
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plot of a SI in Figure [8] for Ma = 1000 GeV, we see that the minimum in a si nearly coincides 
with the zero of the coupling. The minimum (not a zero value) in Figure [8] does not exactly 
coincide with the zero of the coupling, but is shifted closer to Arg{fj,) = 7r/2 because the 
zero value of the real part of the if°-neutralino coupling occurs close to Arg{fi) = 7r/2, as 
shown in Figure iGlFI When (|/i|,Mi) = (175, 110) GeV, the coupling of the lightest Higgs to 
the lightest neutralino has no zero, and asi has no deep minimum, as shown by Figure [51 
For Ma = 200 GeV, the state is much lighter and produces a much larger contribution to 
asi- In this case, the minima of asi are closer to 7r/2, near the zeroes of the if°-neutralino 
coupling, as can be seen in Figure El 

The values of the electron EDM, to be discussed in the next section, are also indicated 
in Figures [3-[9l Among the direct detection experiments, CDMS excludes the region above 
the line labeled as CDMS 2004. The lower lines indicate the projected sensitivities of future 
experiments: CDMS [GSj, ZEPLIN [69] and XENON 



Argiyi) = 




Argm) = %n. 



30 60 70 80 90 100 110 120 130 




50 60 70 SO 90 100 110 120 130 
fGeV) 



Figure 10: Spin independent neutralino-proton elastic scattering cross sections as a function 
of the neutralino mass for Arg{n) = (left) and Arg{n) = 7r/2 (right). Red (dark gray), 
green (medium gray) and yellow (light gray) dots represent models in which the neutralino 
density is above, consistent or below the 2 a WMAP bounds. Hatching indicates the region 
excluded by chargino searches at LEP. The top (blue) solid line represents the 2004 exclusion 
limit by CDMS. The lower solid (cyan) lines indicate the projected sensitivity of CDMS, 
ZEPLIN and XENON, respectively. 

In Figure [TOl we examine the dependence of the direct dark matter detection on the 
phase of /x. In order to do this, we conducted a random scan over the following range of 

^ If the heavy Higgs state is decoupled completely, we find that the minimum of the scattering cross 
section coincides exactly with the zero of the /I'^-neutralino coupling. 
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MSSM parameters: 

-(80 GeV)^ < < 0, 100 GeV < < 500 GeV, 50 GeV < Mi < 150 GeV, 
200 GeV < < 1000 GeV, 5 < tan/? < 10. (17) 

The parameters which are not scanned over are fixed as in Section 2. The result of the scan, 
projected on the stop mass versus neutrahno mass plane, is shown by Figure [lOl Here we 
plot fcTsi as the function of the lightest neutralino mass, where 

rfncDMhy 0.095 a 0.095 > ncDMh'' , . 

^ \ I if 0.095 < QcDMh^ ^ ' 

accounts for the diminishing fiux of neutralinos with their decreasing density [66] For 
models marked by yellow (light gray) dots the neutralino relic density is below the 2 cr 
WMAP bound, while models represented by green (medium gray) dots comply with WMAP 
within 2 a. Models that are above the WMAP value by more than 2 a are indicated by 
red (dark gray) dots. The area indicated by hatching is excluded by the LEP chargino 
mass limit of 103.5 GeV. The top solid (blue) line represents the 2004 exclusion limit by 
CDMS [HZ]. The lower solid (cyan) lines indicate the projected sensitivity of the CDMS [68], 
ZEPLIN [69] and XENON [58j experiments. 

The structure of this scatter plot is clear by examining Figures [UlSl As shown on these 
plots by the gray direct detection contours, the spin- independent cross section, asi, decreases 
for increasing values of Therefore, the low asi region in Figure [TO] is in one to one 
correspondence with the large region in Figures [THHl For large values of the lightest 
neutralino mass is approximately given by Mi, hence, increasing values of Mi in Figures [T]l3] 
correspond to increasing values of m^_^ in Figure [10] and the same annihilation regions, via 

and A^/H^ resonances, and stop coannihilation regions of Figures [HIS] can be identified 
in a clear way in Figure [TUl The LEP excluded, hatched area of m^^ < 103.5 GeV, preserves 
its hyperbolic shape for m^^ < 85 GeV. 

Presently, the region above the (blue) top solid line is excluded by CDMS. In the near 
future, for Arg{fi) = 0, CDMS will probe part of the region of the parameter space where 
the WMAP dark matter bound is satisfied. In this region, due to their enhanced Higgsino 
components, neutralinos mainly annihilate to gauge bosons or, due to the small mass gap, 
they coannihilate with charginos. The ZEPLIN experiment will start probing the stop- 
neutralino coannihilation region together with the annihilation region enhanced by s-channel 

resonances. Finally, XENON will cover most of the relevant parameter space. Prospects 
for direct detection of dark matter tend to be worse for large values of the phase of /x, 
Arg{^) ~ 7r/2. As seen from Figures [7}l9l this phase can lead to cancellations which suppress 
the direct detection cross section. In the event of such a cancellation, a detector with the 
sensitivity of ZEPLIN is needed to start probing the parameter space, and not even XENON 
will be capable of fully exploring this model. 

The experimental limits for dark matter detection rely on the standard assumptions of a dark matter 
flux incident on the earth, based on the observational evidence that points to a roughly spherical distribution 
of dark matter distribution in the galaxy, and a local dark matter velocity comparable to the speed of the 
sun within the galaxy. 
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4 Constraints on CP Violating Phases 



4.1 Electron EDM Constraints 

The MSSM can accommodate many CP violating phases in addition to the CKM phase 
present in the SM. Such phases, however, are very highly constrained by the experimental 
limits on the electric dipole moments (EDM) of the electron, neutron, and ^^^Hg atom. Of 
these, we will focus our attention on the electron EDM since it is the best measured, the 
least plagued by theoretical uncertainties, and for the phases relevant to the model under 
study gives the strongest constraint. The upper bound on the electron EDM comes from 
measurements of the EDM of the "^^^Tl atom. For the phases considered in this work and in 
the absence of Higgs mixing, the CP-odd electron-neutron operator studied in [35] vanishes, 
and the ^os^l EDM is due almost entirely to the electron EDM. This translates into a limit 
on the electron EDM of [72] 

141 < 1-6 X 10"^^ e cm, (19) 

at 90% CL. 

In the MSSM, the leading order contributions to the electron EDM come from one-loop 
diagrams containing an intermediate selectron or sneutrino. For 0{1) phases, these loops 
generate an EDM well above the experimental limit unless these sfermions are taken to be 
quite heavy, > 10 TeV [73]. The neutron and ^^^Hg EDM constraints require that the 
other first and second generation sfermions be very heavy as well. This feature arises in 
several models considered in the literature [ZH USi I76l ED]- Such large first and second 
generation sfermion masses present no problem for EWBG since they couple very weakly 
to the Higgs bosons, and have only a minor effect on the final CP asymmetry ^2]. With 
respect to EWBG, a much more dangerous contribution arises at two-loops. 

At the two-loop order there are relevant contributions to the electron EDM from loops 
containing intermediate charginos and Higgs bosons. Since EWBG demands that the charginos 
be fairly light, < 500 GeV, these contributions cannot be suppressed by taking large 
chargino masses. On the other hand, these terms can be reduced by taking large Ma or 
small tan f3. The phase associated with this contribution comes primarily from the chargino 
mass matrix, which is the same phase that generates the baryon asymmetry, and lower values 
of Ma can enhance the baryon asymmetry. Consequently, the electron EDM bound presents 
a particularly severe constraint on EWBG within the MSSM. 

We have examined whether it is possible for EWBG to generate the observed baryon 
asymmetry while obeying the electron EDM bounds. The two-loop contributions to the 
electron EDM due to intermediate charginos and Higgs were calculated following [Ml [35] . 
The method of [33] was used to calculate the baryon asymmetry generated by EWBG. In 
our analysis, we have fixed M2 = 200 GeV, and varied /x, Arg{fi), tan/3, and Ma- We also 
assume a bubble wall velocity of = 0.05 and a wall width of Lyj = 20/T. Both of these 
values are fairly typical, and tend to maximize the baryon asymmetry generated in the phase 
transition. 

The dependence of the baryon asymmetry (relative to the value needed for big-bang nu- 
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cleosynthesis (BBN)), ti/tibbn, on and Ma is illustrated in Figure [TTl In this plot, we 
have taken the phase to be maximal, sm{Arg{^)) = 1, and have set tan/? = 5. For other val- 
ues of these parameters, the baryon asymmetry scales with sin{Arg{fi)) and (approximately) 
with sin 2(3. There are two main contributions from the CP violating currents of charginos 
and neutralinos to the baryon asymmetry in the MSSM. The first is proportional to the 
change in (3 going from the symmetric phase to the broken phase and exhibits a resonance 
at M2 = |yu|, but is highly suppressed for large values of Ma- The second contribution is 
independent of Ma, and falls off smoothly as becomes large. Both contributions go to 
zero as M2 becomes large. 




Figure 11: Baryon asymmetry generated by EWBG relative to that required by big-bang 
nucleosynthesis for M2 = 200, tan (3 = 5, and sm{Arg{fi)) = 1. 

Figures [T^ a) andfT^b) show the regions in the M4 and My^— tan/3 planes consistent 
with both EWBG and the experimental bound on the electron EDM. Here, we have scanned 
over the ranges 

3 < tan /? < 10, 100 GeV < Ma < 1000 GeV, 100 GeV < < 1000 GeV, (20) 

with M2 = 200 GeV and the rest of the parameters as in Section 2. In Figure [T2k) we see 
that in the allowed region, is confined to the range 110 < < 550 GeV, while M4 must 
be greater than about 200 GeV. The limits on are due to the effect of this parameter on 
the chargino mass. For < 110 GeV, the lighter chargino has mass below the experimental 
bound, m^-^ > 103.5 GeV |17j, while for large EWBG becomes less efficient. The lower 
bound on Ma arises for two reasons. For small Ma the two-loop contribution to the electron 
EDM is enhanced. At the same time the mass of the lightest Higgs is suppressed. The effect 
of the Higgs mass constraint can also be seen in Figure [T^b), in which this bound results in 
a lower limit on tan (3. The allowed region is cut off for larger values of tan (3 since this tends 
to enhance the two-loop contributions to the electron EDM. 
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Figure 12: Parameter regions consistent with EWBG and the electron EDM hmit. In these 
plots, we have taken M2 = 200 GeV and varied Arg{^) over the interval [0, vr]. 

From Figure [121 see that it is possible to generate the baryon asymmetry via EWBG 
in the MSSM while satisfying the experimental constraints on the electron EDM and the 
mass of the lightest Higgs boson. Although this is reassuring, the EWBG scenario is still 
very strongly constrained by the electron EDM. This can be seen in Figure [T3| which shows 
the range of values of dg obtained in our scan that are consistent with EWBG, the current 
electron EDM bound, and the Higgs mass limit. For < 1000 GeV, an order of magnitude 
improvement of the electron EDM bound, \de\ < 0.2 x 10"^''' ecm, will be sufficient to test 
this baryogenesis mechanism within the MSSM. However, we should also point out that 
the calculation of the baryon asjTiimetry from EWBG has 0{1) uncertainties associated 
with the values of the bubble parameters, the wall velocity, and the derivative expansion 
used to derive the diffusion equations. Hence, the limits on EWBG presented here may be 
somewhat more (or less) severe than they really are. Furthermore, we have not considered 
the possibility of fortuitous cancellations between different EDM contributions, for instance 
between the one-loop and two-loop terms (for lighter sfermions), which could further reduce 
the value of the electron EDM. 

4.2 Constraints from BR(6 57) 

The presence of a light stop, light charginos, and a light charged Higgs boson may induce 
relevant effects on flavour changing neutral currents associated with the bottom quark |78]. 
One of the most sensitive experimental measurements of such effects is the branching ratio 
of the decay of a bottom quark into a strange quark and a photon [79]-[82j. A realistic 
calculation of these effects, however, cannot be performed without knowledge of the flavour 
sector of the theory. Even for the large values of the bottom squark masses we consider in 
this work, of order of a few TeV, the contributions coming from the interchange of gluinos 
and down squarks may be as large as the ones coming from the stop-chargino loops [83]. 
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Figure 13: Range of values of the electron EDM for parameter sets consistent with EWBG. 



In the following, we shall present the results for the branching ratio of this rare decay, 
assuming that the only relevant contributions beyond the SM ones are those associated 
with the charged Higgs and stop-chargino loops. While the former tend to increase the 
BR{b 57) compared to the SM value, the latter has a non-trivial dependence on the CP 
violating phase. The experimental value of BR{b 57) is given by [77] . 

BR(6 ^ S7) = (3.54+°i°) x 10'^ (21) 
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Figure 14: BR(6 S7) as a function of Arg{^) for values of the CP-odd Higgs mass 
Ma = 200 GeV (left), and M4 = 1000 GeV (right). The stop parameters were chosen as in 
section 2, and the chargino and neutralino mass parameters are given by (|/i|,Mi) = (350, 
110) GeV (sohd lines), (175, 110) GeV (dashed lines), and (300, 60) GeV (dotted lines). The 
dot-dashed bands represent the present experimental range at the 2 a level. 

In figure [H] we display the value of BR(6 — > 57) as a function of the phase of the Higgsino 
mass parameter, fi, for = 200, 1000 GeV. The stop sector parameters have been chosen as 
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in section 2, while the chargino and neutrahno mass parameters are taken to be Mi) = 
(300, 60) GeV (sohd hnes), (350,110) GeV (dashed hnes), and (175,110) GeV (dotted hnes). 

As is apparent from the figure, in the absence of other sources of fiavour violation, a 
light CP-odd Higgs scalar with mass of about 200 GeV is highly restricted by BR[h —>■ sj). 
Negative values of fiXt, where Xt = At — /i*/ tan (3, are necessary to keep the predicted 
branching ratio close to the experimentally allowed rangej^ This is due to a cancellation 
between the charged Higgs and the squark-chargino contributions to the branching ratio 
when nXf is negative. Otherwise these contributions interfere constructively with each 
other and with the SM contribution. For both signs of {iXt, the branching ratio is largest 
when Arg{ix) = it and smallest for Arg{fi) = 0. Since the branching ratio tends to be 
somewhat high for Ma ~ 200 GeV, even with /xXj < 0, small values of Arg{fi) are preferred 
in this case. 

Larger values of the CP-odd Higgs mass are consistent with the measured value of 
BR(6 — i> 57) over a wide range of values of Mi, fi, and Arg{jj). For moderately large 
values. Ma < 1000 GeV, negative /xXj < is preferred. For Ma > 1000 GeV, the charged 
Higgs contribution decouples leaving only the stop-chargino corrections. These corrections 
tend to give a branching ratio that is near the upper part {^Xt > 0) or lower part [fiXf < 0) 
of the experimentally allowed range for = 700 GeV, as we have considered here. Thus, 
smaller Arg{^) is preferred for fiXt > 0, while Arg{ii) ~ tt is preferred for /iX^ < 0. The 
chargino corrections can be reduced in size by taking slightly smaller values of \Xt\, or by 
invoking small fiavour violation effects in the down squark sector. 

5 Conclusions 

Electroweak baryogenesis provides a mechanism for the generation of the baryon asymmetry 
that relies only on physics at the weak scale. It is therefore testable at high energy physics 
facilities in the near future. In a previous work, we showed that satisfactory dark matter 
abundance may be obtained in the presence of a light stop like the one consistent with 
electroweak baryogenesis, and analyzed the impact of the allowed parameter space for stop 
searches at hadron colliders. No CP violating effects were considered. 

In this work, we have analyzed the effect of CP violating phases, as required for EWBG, 
in conjunction with a light stop, with mass below the top quark mass, and a light Higgs with 
mass below 120 GeV. We have shown that these phases have only a minor impact on the 
stop-neutralino parameter space leading to a consistent relic density. Large phases, however, 
have a relevant impact on direct dark matter detection rates and induce large corrections to 
the electron electric dipole moment. 

We have also shown that, for the phases necessary to obtain an acceptable baryon asym- 
metry and in the limit of heavy squarks, of order a few TeV, the predicted values of the 
electron electric dipole moment tend to lie within an order of magnitude below the reach 

^Recall that if the phases originate from a common gaugino phase and a U(l)ij transformation is used to 
transfer this phase to /i and At, the product pi At remains real but can have either sign. See Eq. dlip . 
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of the present bounds. Even in the case of very heavy squarks, two-loop induced EDM's 
become relevant. Assuming no cancellations between one- and two- loop corrections, one can 
obtain strong bounds on the allowed parameter space: While small values of tan/3 are ex- 
cluded since they lead to unacceptably small values of the Higgs mass, large values of tan P 
tend to lead to unacceptably large values of the electron EDM or small values of the baryon 
asymmetry. On the other hand, for moderate values of tan /5 ~ 7, the Higgs boson mass 
may be large enough to evade the LEP bounds, even for values of as small as 200 GeV. 
For this particular value, and for ~ M2, the baryon asymmetry may be large enough to 
be consistent with observations, even for small values of the phases, of order 0.1, for which 
the EDM's are consistent with the present experimental bounds. 

In the above, we have not discussed the prospects of stop searches at hadron and lepton 
colliders. As discussed in Ref. [36], stop searches become very challenging in the region 
where stop-neutralino coannihilation becomes relevant, both at the LHC and the Tevatron 
collider |84j, due to the small mass difference between the stop and the neutralino. An 
acceptable dark matter density may be obtained for mass differences as small as 20 GeV, for 
which the charm particles proceeding from the stop decay are soft, making the stop detection 
difficult. As shown in Figures 1-3, the presence of CP violating phases doesn't affect this 
result. 

The linear collider signatures of MSSM Baryogenesis have been discussed in Ref. |78j . A 
linear collider represents the best possibility for confirming this scenario since it provides 
the opportunity of performing precise measurements of the chargino system and hence the 
possibility of observing a non-zero phase of the /i parameter [85]. Precise measurements of 
the stop system also become easier at a linear electron-positron collider [86]. For instance, 
the LEP collider was able to set limits on the stops even for a mass difference with the 
neutralino of about 1 GeV. Preliminary studies of stop searches at the linear collider |87j 
show that a 500 GeV ILC may be able to detect a light stop for mass differences as small as 
a few GeV. As said before, in the region of parameters where stop-neutralino coannihilation 
leads to a value of the relic density consistent with experimental results, the stop-neutralino 
mass difference is never much smaller than 20 GeV, and hence an ILC will be able to explore 
this region efficiently 

In summary, the requirement of a consistent generation of baryonic and dark matter in 
the MSSM leads to a well-defined scenario, where, apart from a light stop and a light Higgs 
boson, one has light neutralinos and charginos, sizeable CP violating phases, and moderate 
values of 5 ^ tan f3 < 10. All these properties will be tested by the Tevatron, the LHC and 
a prospective ILC, as well as through direct dark-matter detection experiments in the near 
future. The first tests of this scenario will come from electron EDM measurements, stop 
searches at the Tevatron and Higgs searches at the LHC within the next few years. 
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